Cavity-free Photon Blockade Induced by Many-body Bound States 
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The manipulation of individual, mobile quanta is a key goal of quantum communication; to 
achieve this, nonlinear phenomena in open systems can play a critical role. We show theoretically 
that a variety of strong quantum nonlinear phenomena occur in a completely open one-dimensional 
waveguide coupled to an N-type four-level system. We focus on photon blockade and the creation 
of single photon states in the absence of a cavity. Many-body bound states appear due to the strong 
photon-photon correlation mediated by the four-level system. These bound states cause photon 
blockade, which can generate a sub-Poissonian single-photon source. 

PACS numbers: 42.50.Ct,42.50.Gy,42.79.Gn 



The exchange and control of mobile qubits of infor- 
mation is a key part of both quantum communication 
and quantum information processing. A "quantum net- 
work" is an emerging paradigm [l|, [2] combining these 
two areas: local quantum nodes of computing or end- 
users linked together by conduits of flying qubits. The 
determinisitic approach to the interaction between the 
local nodes and the conduits relies on cavities to provide 
the necessary strong coupling. Indeed, strong coupling 
between light and matter has been demonstrated using 
cavities in both the classic cavity quantum electrodynam- 
ics (QED) systems [3| and the more recent circuit-QED 
implementations [J] . This has enabled the observation of 
nonlinear optical phenomena at the single photon level, 
such as electromagnetically induced transparency (EIT) 
[5|, |6| and photon blockade jTl, Isf . Experiments have also 
demonstrated the efficient exchange of information be- 
tween a stationary qubit (atom) and flying qubits (pho- 
tons) [9l|l0|. However, scaling cavity systems to a multi- 
node quantum network is still challenging because of the 
difficulty of connecting cavities and managing losses. 

A new scheme for achieving strong coupling between 
light and atoms (or artificial atoms) has been recently 
proposed based on one-dimensional (ID) waveguides [lll - 
[l3|, dubbed "waveguide QED" [ll [l|l. Tremendous 
experimental progress in achieving strong coupling has 
occurred in a wide variety of such systems: a metal- 
lic nanowire coupled to a quantum dot [17], a diamond 
nanowire coupled to a nitrogen-vacancy center [18J|_, a 
photonic nanowire with an embedded quantum dot |19j, 
and a ID superconducting transmission line coupled to 
a flux qubit [20[. In these systems, "strong coupling" 
means that the majority of the spontaneously emitted 
light is guided into waveguide modes; it is achieved 
through the tight confinement of optical fields in the 
transverse direction. Furthermore, waveguide systems 
are naturally scalable for quantum networking [ij. The 
key physical element introduced by the waveguide QED 
geometry is that the atom couples to a continuum of 
modes. This relaxes the restriction of working with a 
narrow cavity bandwidth; more importantly, interaction 



with a continuum brings in novel many-body effects that 
have no analogue in a cavity. 

In this work, we show that the nonlinear optical phe- 
nomena EIT, photon blockade, and photon-induced tun- 
neling emerge in a waveguide system for parameters [20| 
that are currently accessible. For these dramatic and 
potentially useful nonlinear effects, it is necessary to 
consider a four-level system (4LS) rather than simply a 
two- or three-level system. Photon blockade and photon- 
induced tunneling have a completely different origin here 
from the cavity case ITI. 12111 : they are produced by many- 
body bound states I^llj,ll5|,l23|, whose amplitude decays 
exponentially as a function of the relative coordinates of 
the photons. Such states do not exist in cavities because 
a continuum of modes in momentum space is needed for 
the formation of bound-states in real space. We demon- 
strate the capability of such a system to generate a single- 
photon source, which is crucial for quantum cryptogra- 
phy and distributed quantum networking. Our work thus 
opens a new avenue toward the coherent control of light 
at the single-photon level based on a cavity-free scheme. 

Motivated by recent experimental advances [2(f\, we 
consider an TV- type 4LS [23|, |24| coupled to a continuum 
of modes in a ID waveguide. Figure 1 shows both a 
schematic and a possible realization using superconduct- 
ing circuits [24j . The Hamiltonian of the system is [12, uA \ 
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where a]^ r{x) s-re the creation operators for left, right- 
going photons at position x and c is the group velocity of 
photons. In the 4LS, the energy reference is the energy 
of state |1) (the ground state), and 62 = UJ21, £3 = £2 — A, 
and £4 = £3 + a;43, where UJ21 and UJ43 are the |1) -H> |2) 
and |3) O |4) transition frequencies, respectively. In the 
spirit of the quantum jump picture [25|, an imaginary 




wave-vector k) and the 4LS initially in its ground state, 
the transmitted part of the single-photon eigenstate is 



FIG. 1: Schematic diagram of the waveguide system and EIT. 
(a) A'^-type four-level system interacting with the waveguide 
continuum. The transitions |1) -f-!- |2) and |3) -f-^- |4) are cou- 
pled to the waveguide modes with strength V, detuning S. 
The transition |2) -H- |3) is driven by a semiclassical control 
field with Rabi frequency Q and detuning A. Here, uj and 
(jjc are the frequencies of the incoming photons and the con- 
trol field, (b) Schematic of possible experimental setup based 
on superconducting charge qubits [24|. Two charge qubits, 
each formed from two Josephson junctions [dark gray (green) 
regions] , are coupled capacitively to each other and to a trans- 
mission line (ID waveguide), (c) Single- photon transmission 
probability as a function of incident photon detuning. EIT 
occurs when the control field is on (f2 = 1.6, solid); when it is 
off, the 4LS becomes a two-level system and EIT disappears. 
Here, the effective Purcell factor is P = 6. 



term is included in the 4LS to model the spontaneous 
emission of the excited states at rate Tj to modes other 
than the ID waveguide continuum. We have assumed a 
linear dispersion and a frequency-independent coupling 
strength V for the relevant frequency range [l^. The 
decay rate into the waveguide modes is F = 2V'^ /c from 
Fermi's golden rule. Below, we assume that level |3) is 
metastable (F3 = 0) and levels |2) and |4) have the same 
loss rate (F2 = F4). 

A figure of merit to characterize the coupling strength 
is given by the effective Purcell factor, P = F/F2. In an 
experiment with surface plasmons coupled to a quantum 
dot 17[, P = 1.5 was achieved. In more recent experi- 
ments, even larger Purcell factors, P = 3 and P > 9, were 



demonstrated with a superconducting transmission line 

19|,l26|, respectively. 



[20| and a GaAs photonic nanowire 
These recent dramatic experimental achievements sug- 
gest that the large Purcell-factor physics, which we now 
discuss, is presently within reach experimentally. 

To study interaction effects during photon transmis- 
sion, we obtain an exact solution of the scattering prob- 
lem defined by Eq. ^ . The scattering eigenstates are 
obtained by imposing an open boundary condition and 
requiring that the incident photon state be a free plane 
wave, an approach adopted previously to solve an inter- 
acting resonant-level model |22| and a two-level system 
problem [141] . For an incident photon from the left (with 
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where tk is the transmission coefficient. In the two- 
photon scattering eigenstate, the transmitted wave is 
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where 61.2 and 71^2 (> 0) are functions of system parame- 
ters (see the Supplementary Information) . The first term 
of (j^RR corresponds to transmission of the two photons as 
independent (identical) particles with the momentum of 
each photon conserved individually. The second term is 
a two-body bound state — note the exponential decay in 
the relative coordinate \x2 —xi\ — with two characteristic 
binding strengths 71 and 72. Such a state results from 
the nonlinear interaction between photons mediated by 
the 4LS. Physically, it originates from the momentum 
non-conserved processes of each individual photon (with 
conservation of total momentum). A similar bound state 
has been found in a ID waveguide coupled to a two- level 
system I3, [ij] , a A- type three- level system [la , as well 
as an open interacting resonant level model [22l |. 

We evaluate the transmission and reflection probabili- 
ties using the S-matrices constructed from the exact scat- 
tering eigenstates [iJ] . Because any state containing a fi- 
nite number of photons is, in practice, a wave packet, we 
consider a continuous mode input state [27| , whose spec- 
trum is Gaussian with central frequency loq and width 
a. Throughout this paper, we set the loss rate as the 
reference unit for all other quantities: F2 = r4 = 1. The 
Purcell factor becomes P = F/F2 = F. For all the nu- 
merical results shown, we take the detuning of the control 
field to be zero, A = 0, and choose a = F2/5 = 0.2. In 
addition, we assume that the transitions |1) O |2) and 
|3) O |4) are at the same frequency {uj2i — W43). 

Figure 1 (c) shows that the transmission probability Ti 
of a single photon has a sharp peak as a function of its 
detuning S = ujn— UJ21, demonstrating the familiar EIT 
phenomenon [l5|, |23| produced by interference between 
two scattering pathways. The width of the EIT peak is 
~ f2 for P ^ 1 and is mainly determined by the control 
field strength Q. Here, we use P = 6, a conservative value 
given the recent advances in experiments [l9|, |20|, l26l |. 
When the control field is turned off, the 4LS becomes a 
two-level system, which acts as a reflective mirror |13j . 

The EIT picture changes dramatically when there are 
two or more photons injected into the system. The 4LS 
mediates an effective photon-photon interaction, which 
in turn affects the multi-photon transmission. We de- 
fine T2 (Ts) to be the two- (three-) photon transmission 





FIG. 2: Photon blockade in transmission, (a) Photon block- 
ade strength P21 (solid) and P31 (dashed) as a function of the 
decay rate F into the waveguide modes (photon blockade is 
strongest when P21 is smallest). The photons are on resonance 
with the 4LS {5 — 0) and Q = 1.6. Strong coupling produces 
substantial photon blockade despite EIT being present in the 
single photon properties, (b) P21 and P31 as a function of the 
incident photon detuning 5 for F = 6 and fl = 1.6. There is 
clear photon blockade on resonance as well as photon-induced 
tunneling {P2i,P3i > 1) on the sides of the resonance. 



probability of the two- (three-) photon scattering pro- 
cess. The strength of the photon blockade P21 is given 
by the conditional probability T2/T1 for transmitting a 
second photon given that the first photon has already 
been transmitted, normalized by the single-photon trans- 
mission probability Ti: P21 = T2/T1. For independent 
photons, there is no photon blockade and P21 = 1- In 
the opposite limit of strong photon blockade, P21 is sup- 
pressed towards zero. Similarly, we define P31 = T^/T^ 
to quantify photon blockade in the three photon case. 

A pronounced photon blockade is shown in Fig. 2(a) 
in the strong coupling regime: the single photon EIT ef- 
fect does not carry over to the multi-photon case. For 
increasing coupling strength, both P21 and P31 approach 
zero. Such photon blockade regulates the flow of photons 
in an ordered manner, enabling coherent control over the 
information transfer process in our cavity-free scheme. 
Taking P = 9 achieved in Ref.|l9|, we obtain the values 
P21 ^ 30% and P31 ~ 7%, showing that the effects pre- 
dicted here are already within reach of experiments. 

The photon blockade occurs despite being in the EIT 
regime: as shown in Fig. 2(b), both P21 and P31 are sup- 
pressed within the EIT window, whose width is set by 
the control field strength J7. However, away from the 
EIT window, P21 and P31 become larger than 1, signal- 
ing a new regime of multi-photon transmission — photon- 



FIG. 3: Many-body bound-state effect causing photon block- 
ade, (a) Transmission as a function of the decay rate F into 
the waveguide modes for n — 1.6. The presence of the many- 
body bound states suppresses the multi-photon transmission, 
(b) Transmission as a function of incident photon detuning 5 
for F = 6 and Q, = 1.6. The bound states strikingly oppose 
multi-photon EIT; off-resonance, their constructive interfer- 
ence produces photon-induced tunneling. 



induced tunneling [2lJ. Previously, photon blockade and 
photon-induced tunneling have been observed in cavity- 
QED and circuit-QED systems 0, |1, [2l|, where the un- 
derlying mechanism is the anharmonicity of the spectrum 
caused by the atom-cavity coupling. We emphasize that 
such anharmonicity is absent in the present cavity-free 
scheme: the photon blockade and photon-induced tun- 
neling here must be caused by a different mechanism. 

To understand the origin of these phenomena, we sepa- 
rate the two-photon transmission probability T2 into two 
parts (Fig. 3): (T2)pw = T^ is the contribution from 
independent transmission and (T2)bs is the contribution 
from both the bound-state term in Eq. (|3]) and the inter- 
ference between the plane-wave and bound-state terms. 
{T2)bs is the result of the many-body interactions in the 
waveguide and is absent for cavities. Similarly, r3 can 
be separated into (T3)pw = T^ and (T3)bs- Figure 3(a) 
shows that, when the photons are on resonance, (T2^3)bs 
is always negative, suppressing the overall transmission. 
The cause of the observed photon blockade is, thus, the 
destructive interference between the two transmission 
pathways: passing by the 4LS as independent particles 
or as a composite particle in the form of bound states. 
Within the EIT window, this conclusion always holds (see 
Fig. 3(b)) — both two- and three-photon transmission are 
strongly suppressed by the many-body bound-state ef- 
fect. In contrast, away from the EIT window, (T2^3)bs 
changes sign and becomes positive (see the Supplemen- 
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FIG. 4: Non-classical photon states. Photon number statis- 
tics quantified by logj^Q(P„/P„,Poiss3on), where the probability 
of n photons in the output state is P„ and Pn.Poisson is for 
a coherent state with the same mean photon number. Pan- 
els (a)-(d) show results for n = 0,1,2, and 3, respectively, 
as a function of the control field strength Q. and the decay 
rate F. The dashed line is a guide to the eye indicating equal 
probabilities, Pn/Pn, Poisson = 1- Note the large range of pa- 
rameters for which the single photon probability is enhanced 
while the multi-photon content is suppressed — an improved 
single photon source. 



tary Information for a detailed analysis of the interfer- 
ence). Destructive interference changes to constructive 
interference, producing photon-induced tunneling. 

As an application, we now show that the 4LS can gen- 
erate non-classical photon states. We assume that the 
4LS is in its ground state initially. We consider an in- 
cident continuous-mode coherent state of mean photon 
number n — 1.0, on resonance with the 4LS. The pho- 
ton number statistics in the transmitted field can be ob- 



tained using the S-matrix method IJ] and is presented 
in Fig. 4 by taking the ratio of the photon number dis- 
tribution of the output state Pn {n — 0, 1, 2, 3) to that of 
a coherent state Pn,poisson having the same mean photon 
number. It is remarkable that, in most of the parame- 
ter space, we have Pi > Pi,poisson while P2 < P2,Poisson 
and P3 < PsjPoisson- This gives rise to a sub-Poissonian 
single-photon source: for example, for T = 9 (P = 9) 
and f} = 1.6, we have Pq = 64% and Pi = 34% with the 
multi-photon probability P„>2 less than 3%, in compar- 
ison with P„>2 ~ 26% in the input. This single-photon 
source comes about because, under EIT conditions, a sin- 
gle photon passes through the system with high probabil- 
ity, while multi-photon states experience photon block- 
ade caused by the bound-state effect. A systematic way 
of improving the quality of the single-photon source is to 
let a coherent state with a large mean photon number 
pass through multiple 4LS devices in series with Faraday 
isolators inserted between each stage. 



In summary, we present a cavity-free scheme to real- 
ize a variety of nonlinear quantum optical phenomena — 
including EIT, photon blockade and photon-induced 
tunneling — in a ID waveguide. Photon blockade and 
photon-induced tunneling have a distinctly different ori- 
gin here compared to the cavity case: a many-body 
bound-state effect. Furthermore, we outline how to use 
EIT and photon blockade in this system to produce a 
single-photon source on demand. On the one hand, the 
demonstrated ability to control the flow of light quanta 
using EIT and photon blockade is a critical step towards 
the realization of an open quantum network. On the 
other hand, the strong photon-photon interaction me- 
diated by the 4LS provides a new candidate system to 
study strongly correlated ID systems, one complemen- 
tary to condensed-matter systems. 

The work of HUB and HZ was supported by the 
U.S. OfSce of Naval Research. 
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Hamiltonian 



The Hamiltonian describing the interaction between a ID continuum and a 4LS at the spatial origin is given by 

H = j dkhujka\k)a{k)+HiLS+ j dkhgia){k)[Si2 + S'n\+h.c\, 
Hals = JZ^U-'^)\3){j\ + ^{\'2){-i\+h.c.) (S4) 

where g is the coupling strength, 5*12 = |1)(2| and 5*34 = |3)(4|. Linearizing the dispersion around the transition 
energy of the 4LS and treating the left-going and right-going photons as separate fields, we obtain a two-mode 
Hamiltonian. Transformation to the real space field operators ai(/j)(a;) = {1/\/2t:) ^ dke^^'^ai^(^fj^{k) yields the real- 
space Hamiltonian 

H^ ldxhc[a}j^{x){-i — )aR{x)-a\{x){~i — )aL{x)]+HiLS+ j dxhv[[a)j^{x) + a\{x)\[Si2 + S34\ + h.c.\, (S5) 



where V = \/2T:g. Finally, the coupling between the 4LS and general modes other than the waveguide can be 
incorporated in the same manner as the waveguide-4LS interaction. Such interaction with general reservoir modes, 
such as vacuum modes, leads to intrinsic dissipation. This term can be treated by the Wigner-Weisskopf approximation 
[3| and so replaced by a non-Hermitian damping term in H^ls |2lj|j]' thus giving rise to the Hamiltonian in Eq. (f) 
of the main text. 

Scattering Eigenstates 

A key insight to solve for the scattering eigenstates is gained by performing the operator transformation 

al/,{x)^[alix)±al{-x)]/V2. (S6) 

The Hamiltonian is then decomposed into even and odd modes: H = H^ + Hq with 

/T 
dx{~i)hcdl{x) — ae{x) + H^ls 

+ f dx hV2V 6{x) [al (x) iSi2 + S3i) + h.c.], (S7) 

dx(—i)hcdl(x)-—do(x). 
dx 

Because [H, he + n4Ls] — [H, ho] — for the number operators hf./o = J dxd\,{x)de/o{x) and n4LS = X]i=2 3 b) 01 + 
2|4)(4|, the total number of excitations in both the even and odd spaces are separately conserved. We will now focus 
on finding the non-trivial even-mode solution and then transform back to the left/right representation. A general 
n-excitation state {n = Ue -\- ?t.4ls) in the even space is given by 

|*n> = ndx^g^^\x)dlixi)---dlixn) 

+ Jdx-'[ft\x)St2 + ft\^)Srs] al(xi)---at(:r„_i), 

+ Jdx^-^f!C\x)S+,dl{x,)---dl{xn-2))\0,l} (S8) 



where |0, 1) is the zero photon state with the 4LS in the ground state |1). 

The scattering eigenstates are constructed by imposing the boundary condition that, in the incident region 
Xi , • • • ,Xn < 0, 5*^"^ {xi , • • • ,Xn) is a free plane wave [5|, |6[ . The one-photon scattering eigenstate with eigenenergy 
_E = cfc is given by 
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[cfc - £2 + A + ir3/2] [cfc - £2 + «r2/2 + ir/2] - ^2/4 ' 
where 6{x) is the step function and F = 21/2/^ is the spontaneous emission rate to the ID continuum. Notice that 



the amphtude fl vanishes because it takes two excitation quanta to excite level |4). Therefore, for one-photon 



scattering, the 4LS behaves as a three-level system and Eq. (|S9p coincides with previous studies [7|, |8[. Transforming 
from the even/odd back to the left/right representation gives rise to the wavefunction of one photon transmission in 
the single-photon scattering, shown in Eq. (2) in the main text, where t^ — [tt + l)/2. 

For two-photon scattering, starting from a free plane wave in the region a;i, X2 < 0, we use the Schrodinger equation 
implied by H^, in Eq. (|S7[) to find the wave function first in the region xi < < a;2 and then for < a;i,X2. The 
method is explained in detail in the appendix of Ref. [6| for scattering from a two level system. We obtain the 
two-photon scattering eigenstate with eigenenergy E = cfci -I- cfc2 
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FIG. 5: The integrand function TB{k\, k2) of the interference term as a function of fciand k2. (a) 5 = 0.0. (b) 5 — 2.0. 

Again, performing an even/odd to left/right transformation [6| , we obtain the wavefunction of two photon transmission 
in two-photon scattering, shown in Eq. (3) in the main text. There, 71/2 = — Re[Ai/2] and bi/2 oc Ci/2- In the case 
that A = Ta == 0, C = |x|6l(-x)/2 and rj = \x\0{x)/2 with x = ^^ - (r + r2)^ /4, where again 0{x) is the step function. 
From the scattering eigenstates, scattering matrices can be constructed using the Lippmann-Schwinger formahsm 
[9|. The output states are then obtained by applying the scattering matrices on the incident photon states (single-, two- 
, and three-photon number states as well as coherent states). By performing measurement operations on the output 
states, we obtain the transmission and reflection probabilities and photon number statistics. A general procedure is 
outlined in Ref. [yj. 

Phase Analysis of Photon Blockade and Photon-Induced Tunneling 

In Fig. 3 of the main text, we show that the observed photon blockade and photon-induced tunneling orignate 
from destructive and constructive interference effects, respectively. In this section, we present a detailed analysis of 
the relative phase between the plane-wave and bound-state terms. The two-photon transmission probability T2 for 
the two-photon scattering is given by 



r2- f dkidk2\T{ki,k2) + B{ki,k2)\^ 



(S12) 



where fci and fc2 are the momenta of the two photons in the wavepacket. The plane-wave term T(fci,/c2) and the 
bound-state term B{ki, k2) take the form 



T{ki,k2) = a{ki)a{k2)tkfy. 



B{ki,k2) - y ^( , \ + , \ \ f dk'a{k')a{ki+k2-k')Cj{k',ki + k2-k'), 
■tl^2 ^ V^i-*7j k2-iljJJ 



(S13) 



where a{k) is the amplitude of a Gaussian wavepacket with central frequency loq and width a. In momentum space, 
the central momentum is ko — uiq/c and the width is a/c, where c is the group velocity of photons in the waveguide. 
By defining the phase difference 9{ki,k2) between T{ki,k2) and B{ki,k2), T2 can be written as 



T2 = f dkidk2\T{ki,k2)\^ + f dkidk22\T{ki,k2)B{ki,k2)\cos9{ki,k2) + f dkidk2\B{ki,k2)\^, (S14) 
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where the first term is the plane-wave term, the second term is the interference between the plane-wave and bound- 
state terms, and the third term is the contribution from the bound-state term. In the main text, the first term is 
denoted {T2)pw and the second and third terms together are called {T2)bs- Denote the intergrand function of the 
interference term as TB{ki, ^2), 

TB{ki, fca) = 2|T(fci, fc2)B(fci, fc2)|cos0(fci, fca) (S15) 

The phase d{ki,k2), or specifically, the sign of cos6'(fci, ^2), determines whether the interference is constructive or 
destructive. 

We numerically evaluate TB{ki,k2) in two cases: 5 — loq ~ W21 — 0.0 and 5 — ujq ~ UJ21 — 2.0. Here, the unit 
of detuning S is set by the loss rate r2. All the other system parameters are the same as in Fig. 3b in the main 
text. Fig. [5] shows TB{ki,k2) as a function of fci — fco and ^2 — fco- As expected for a Gaussian packet, the value 
of TB{ki,k2) is centered at ki = k2 — ko in both cases. However, the sign of the peaks in the two cases differs. 
For S — 0.0, a negative peak indicates destructive interference, giving rise to photon blockade when the incident 
photons are on resonane with the 4LS. For S = 2.0, a positive peak indicates constructive interference, producing 
photon-induced tunneling when the photons are far off resonance. 
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